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Avoidance of future singularities in loop quantum cosmology 



O 

o 

(N 



M. Sami,^'^ Parampreet Singh,"'* and Shinji Tsujikawa^ 

Centre for Theoretical Physics, Jamia Millia, New Delhi-1 10025, India 
^Department of Physics, Jamia Millia, New Delhi- 110025, India* 
'^Institute for Gravitational Physics and Geometry, 
Physics Department, Penn State, University Park, PA 16802, U.S. A} 
'''Department of Physics, Gunma National College of Technology, Gunma 371-8530, Japan} 

(Dated: February 6, 2008) 

We consider the fate of future singularities in the effective dynamics of loop quantum cosmology. 
Non-perturbative quantum geometric effects which lead to p'^ modification of the Friedmann equa- 
tion at high energies result in generic resolution of singularities whenever energy density p diverges 
at future singularities of Friedmann dynamics. Such quantum effects lead to the avoidance of a Big 
Rip, which is followed by a recoUapsing universe stable against perturbations. Resolution of sudden 
singularity, the case when pressure diverges but energy density approaches a finite value depends 
on the ratio of the latter to a critical energy density of the order of Planck. If the value of this ratio 
is greater than unity, the universe escapes the sudden future singularity and becomes oscillatory. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



Observations suggest that the current universe is dom- 
inated by a matter component which leads to an acceler- 
ated expansion of the universe- called dark energy (see 
Refs. [1] for review). This has stimulated a study that 
our universe may face a future singularity. Such future 
singularities typically arise if the universe is dominated 
by matter which violates dominant energy condition and 
causes a state of super-acceleration of the universe be- 
fore leading it to a singularity. They can occur due to 
divergence either of the energy density p and/or the pres- 
sure density p of the matter content. For example if the 
universe is filled with a phantom dark energy with a con- 
stant equation of state w less than — 1 [2] , this leads to a 
Big Rip singularity at which both p and p diverge with 
a finite time [3] . Barrow pointed out a possibility to ob- 
tain a sudden future singularity at which p is finite but 
p diverges [4]. Depending on the equation of state of 
dark energy, future singularities have been categorized in 
diff'erent classes [5]. 

Existence of future singularities in Friedmann- 
Robertson- Walker (FRW) cosmology reflects the vulner- 
ability of standard Friedmann dynamics whenever p or p 
become of the order of Planck values. This indicates that 
limit of validity of general relativity has been reached and 
inputs from quantum gravity are necessary to probe the 
dynamics near the singularity. Resolution of singularities 
using Wheeler-DeWitt quantization has been attempted 
[6] but has met with little success. One of the primary 
reasons for its failure has been a lack of a fundamen- 
tal theory which can guide quantization in the Wheeler- 
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DeWitt framework. Issue of resolution of past [7] and 
future singularities [8] have been investigated using per- 
turbative corrections in string theoretic models. These 
analysis indicate that generic resolution of singularities 
may only be accomplished using non-perturbative correc- 
tions. In particular in the absence of an analysis which 
uses non-perturbative quantum gravitational modifica- 
tions to model the dynamics of dark energy, the fate of 
future singularities has remained an open problem. 

Loop quantum gravity (LQG) is a leading background 
independent non-perturbative quantization of gravity [9] 
which has been very well understood in the cosmologi- 
cal setting in loop quantum cosmology (LQC) [10]. To 
its success, LQG has dealt with various singularities in 
cosmological setting [11-15] and techniques have also 
been used to resolve singularities in black hole space- 
times [16]. Recent investigations have revealed that non- 
perturbative loop quantum effects lead to a p^ modifi- 
cation of the Friedmann equation with a negative sign 
[13; 17, 18]. The modification becomes important when 
energy density of the universe becomes of the same order 
of a critical density pc. The resulting dynamics generi- 
cally leads to a bounce when our fiat expanding universe 
is evolved backwards [12-15]. 

Since important insights have been gained on resolu- 
tion of space- like singularities in LQC, it offers a natu- 
ral arena to investigate the fate of future singularities. 
This is the goal of the present work. Using the effective 
Friedmann dynamics which has emerged from LQC we 
would analyze the way non-perturbative quantum gravi- 
tational effects modify the dynamics near future singular- 
ities. The plan of this paper is as follows. In the next sec- 
tion we would briefiy review the way an effective modified 
Friedmann dynamics is obtained from the discrete quan- 
tum dynamics in LQC. Since, Wheeler-DeWitt quantiza- 
tion, which like LQC is a mini-superspace approach, has 
been unsuccessful in resolving space-like singularities we 
would highlight some differences which emerge with LQC 
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(for details see Ref. [13]). In Sec. Ill we analyze in detail 
the fate of three types of future singularities- type I (the 
Big Rip): scale factor a, energy density p and pressure 
p becoming infinite in finite time, type II (sudden): p 
becoming infinite with finite p in finite time, and type 
III: p and p diverging with finite a in finite time. We 
will show the LQC can successfully resolve type I and 
type III singularities for generic choice of initial condi- 
tions. Resolution of type II singularities though depends 
on the amplitude of model parameters. We conclude with 
a summary of our results in Sec. IV. 



II. EFFECTIVE DYNAMICS IN LOOP 
QUANTUM COSMOLOGY 

In LQG the phase space of classical general relativ- 
ity is expressed in terms of SU(2) connection A\ and 
densitized triads Ef. In loop quantum cosmology (LQC) 
[10], due to underlying symmetries of the FRW spacctime 
the phase space structure simplifies and can be castcd in 
terms of canonically conjugate connection c and triad p 
which satisfy {c, p} = k7/3, where k = SttG (G is grav- 
itational constant) and 7 is the dimensionless Barbero- 
Immirzi parameter (which is set by the black hole ther- 
modynamics in LQG, as 7 w 0.2375). On the space of 
physical solutions of general relativity they are related to 
scale factor and its time derivative as: c = 7a and p — a^. 

The elementary variables used for quantization in LQC 
are the triads and holonomies of connection over edges 
of loops: hi{p) = cos(/^c/2) 4- 2 sin(/ic/2)ri, where are 
related to Pauli spin matrices as Ti = —ioijl and p is 
related to the length of the edge over which holonomy is 
evaluated. The algebra generated by holonomies is that 
of almost periodic functions of c with elements of the 
form: exp(z^c/2). On quantization, though holonomies 
have well defined quantum operators there are no quan- 
tum operators for c in LQC (as in LQG). The kinematical 
Hilbert space in LQC is = -^^(MBohr, ^A^Bohr) where 
^Bohr is the Bohr compactification of the real line and 
/^Bohr is the Haar measure on it. Note that the Hilbert 
space is different from the one in Wheeler-DeWitt quan- 
tization: TiwDW = L'^{^,dp). The triad and thus the 
scale factor operator in LQC have a discrete eigenvalue 
spectrum and quantum constraint, obtained by express- 
ing the classical constraint in terms of holonomies and 
positive powers of triad and then quantized, in LQC leads 
to a discrete quantum difference equation whose all so- 
lutions are non-singular - another important distinction 
from the Wheeler-DeWitt theory. 

Physical predictions can be extracted from LQC by 
construction of a physical Hilbert space. By identifying 
Dirac observables on this space, information about dy- 
namics can be extracted using ideas of emergent time. 
On constructing coherent states we can then find out the 
expectation values of Dirac observables and compare the 
quantum dynamics with the classical one. It turns out 
that when a flat expanding universe is evolved backward 



using loop quantum dynamics, instead of ending in big 
bang singularity it bounces at Planck scale to a contract- 
ing branch [12-14]. 

The coherent states used to analyze the details of quan- 
tum dynamics also play an important role in obtaining an 
effective Hamiltonian description of dynamics governed 
by quantum difference equation. This can be done by 
using methods of geometric formulation of quantum me- 
chanics [19] where one notes that quantum Hilbert space 
can be regarded as a quantum phase space with a bun- 
dle structure. The classical phase space forms the base 
of this bundle, whereas fibers consist of states with same 
expectation values of conjugate variables. Horizontal sec- 
tions of the bundle are isomorphic to the classical phase 
space. Using coherent states one can then find horizontal 
sections which are preserved by quantum evolution which 
then leads us to an effective Hamiltonian with loop quan- 
tum modifications [20, 21]: 
3 

Coff = TTTj asin^(/2c) -I- Cm . (1) 
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Here p is the kinematical length of the edge of a square 
loop which has the area given by the minimum eigenvalue 
of the area operator in LQG [18] and Cm corresponds 
to matter Hamiltonian which in general contains modifi- 
cations due to rcgularization of inverse scale factor [10]. 
These modifications are negligible for large universes and 
would not be considered in the present work. 

The modified Fricdmann equation can then be ob- 
tained by using the Hamilton's equation: 

P = iP,Ccff| = — y^) 
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in the effective Hamiltonian constraint Cofr ~ 



Pc 



(3) 



with Pc = V^/{167r-/^G^h) [18] where h is Planck con- 
stant. Along with the conservation law: 



p + 3H{p + p)=0, 



(4) 



Eq. (3) provides an effective description of Fricdmann dy- 
namics which very well approximates the underlying dis- 
crete quantum dynamics and confirms with the picture of 
bounce which occurs when p = pc (of the order of Planck 
density) [13-15]. In the classical limit S — *■ one has 
Pc 00, thus classically non-singular bounce is absent. 
Further, for p Pc the modified Fricdmann equation 
reduces to the standard one. Interestingly p^ modifica- 
tions also appear in string inspired braneworld scenarios 
and it turns out that there exist interesting dualities be- 
tween two frameworks [17]. However, such modifications 
in braneworlds usually appear with a positive sign and 
a bounce is absent (unless one assumes existence of two 
time- like dimensions [22] ) . We will now address the issue 
of future singularities in the effective dynamics of LQC. 
Our treatment of matter which leads to such singularities 
would be phenomcnological and at an effective level (as 
in Ref. [23]). 
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III. AVOIDANCE OF FUTURE 
SINGULARITIES 

For the analysis of the fate of future singularities in 
LQC, it is useful to first obtain the rate of change of 
Hubble parameter from Eqs. (3) and (4): 



2 \ Pc 



(5) 



where w is the equation of state: w = p/ p which in 
general may not be a constant. 

It is convenient to define two variables 



y 



P 



~ 3iJ2 
Then from Eq. (3) we find 

y = l- l/x. 



(6) 



(7) 



Since is positive, the variables y and x are in the 
ranges < y < 1 and x > 1. From Eqs. (4) and (5) we 
obtain the differential equation for the variable x: 



dx 
d/V 



-'d{l + w)x{x - I) , 



(8) 



where = In (a). 

Since H can change the sign, it will be convenient to 
solve differential equations in terms of a cosmic time t 
rather than N . Defining two dimensionless quantities 
t = Hct and H = H/Hc, where He = \/ k/9c/3, we find 
that Eqs. (4) and (5) are written as 



at 



-(l + «;)y(l-2y). 



AH 

"dT ~ ^ 2~ 
together with the constraint equation 

H'' = v{\-y). 

Combining equations (9) and (10) gives 



(9) 
(10) 

(11) 



|| = ^(l + HV(3-4y)-3il2;^. (12) 

We will study several equations of state which, in stan- 
dard Einstein gravity, give rise to various types of fu- 
ture singularities [5] . Our interest is to clarify the role of 
loop quantum modifications on the following singularities 
which arc known to exist in standard Einstein gravity: 

• Type I ("Big Rip") : For t ^ tg, p ^ oo, |p| — > oo, 

H oo and a oo 

• Type II ("sudden") : For t ^ tg p ^ ps, \p\ oo, 

H —f Hg and a — > 

• Type III : For t tg, p oo, \p\ — > oo, — > oo 
and a — > Cs . 



Here tg, Ps, Hg and ag are constants. The type I singu- 
larity appears for constant w less than —1 [3]. The type 
II is a sudden future singularity [4] at which p and a are 
finite but p diverges. The type HI appears for the model 
with p = — p — Ap^^ with (3 > 1 [24]. In what follows we 
shall study each case separately. 



A. Type I singularity 

Let us consider a constant equation of state, w. In this 
case Eq. (8) is easily integrated to give 



1 



1 - Aa-3(i+'") 



y 



^a-3(i+"') , (13) 



where A is a positive constant. 

When w > —I the solutions in an expanding uni- 
verse approach the fixed point {x,y) = (1,0), which cor- 
responds to the standard Einstein gravity. Meanwhile 
when ui < — 1 one has x — > oo and y 1 with a 
scale factor satisfying Aac^'"^^^'' = 1. In this case p 
approaches a constant value pc b,s a —^ Uc- From Eq. (3) 
we find that the Hubble parameter becomes zero at this 
point. This equation [or equivalently Eq. (11)] also tells 
us that the Hubble parameter H varies between its max- 
imum and minimum values given by H^^j^ = 1/2 and 
ffmin = ^1/2 respectively. We also notice from Eq. (10) 
that dH/dt < at the time when y becomes y = 1, 
leading to the decrease of H. Equation (9) tell us that 
y begins to decrease after it has reached its maximum 
value y = 1 corresponding to H = 0. 

We can now understand the fate of the universe for 
u" < — 1 in LQC. A qualitative description of the evolu- 
tion can be obtained by using Eqs. (9), (10) and (11). 
Let us begin to examine the evolution with a positive 
initial value of H. From Eqs. (9) and (10), we find that 
both y and H grow until y reaches y = 1/2. Such a 
behavior of p and H is generic to a phantom dominated 
universe with constant w in the standard FRW cosmol- 
ogy which then leads to a Big Rip. The LQC correction 
changes this cosmic evolution in a crucial manner. The 
Hubble rate begins to decrease after it reaches a maxi- 
mum value at y = 1/2, whereas y continues to grow until 
H drops below zero [see Eqs. (9) and (10)]. As explained 
above, the variable y starts to decrease with dH/dt < 
after it has reached its maximum value y = 1. When y 
becomes smaller than 1/2, dH/di changes its sign after 
which _ff (< 0) increases toward 0. This stage corresponds 
to a rccoUapsing universe that asymptotically approaches 
the fixed point {x,y) = (1,0). From Eq. (12) together 
with Eq. (11) we find that the asymptotic behavior is 
given by y oc t^^ and H oc — t^^. Note that for negative 
H the fixed point {x,y) = (1,0) is stable against per- 
turbations as can be checked by linearly perturbing the 
system (9) and (10). In fact two eigenvalues of a matrix 
for perturbations [25] are and —3(1 + iv)H, where the 
latter is negative ior H < 0. 
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FIG. 1: Evolution of the Hubble parameter and the variable 
y = P/Pc for w — —1.5 with initial conditions j/i = 0.01 and 
Hi = [yi{l — y,;)]^''^. The Big Rip singularity is avoided in 
the presence of loop quantum modifications to the Friedmann 
dynamics. 

We have numerically solved Eqs. (9) and (10) for 
If < — 1 with initial conditions H > 0. In Fig. 1 we plot 
an example for the evolution of H and y when w = —1.5. 
Our numerical results clearly confirm the qualitative be- 
havior of the evolution presented above. 

Thus we have shown the Big Rip singularity is beau- 
tifully avoided in the frame work of LQC. The solutions 
finally approach a contracting universe in standard Ein- 
stein gravity {H and p-^Oasi— >oo). Wc note 
that when w > —1 bouncing solutions can be obtained if 
H <0 initially [12-15]. 

B. Type II singularity 



where B = and r = po/Pc- Integrating this 

equation gives 

(^r-l + i^ =(^r-l + ^^ -B(7+l)iV, (16) 

where we chose the initial condition x = Xi at N — 0. 
This shows that x gets larger with the increase of N, in 
which case y = p/pc grows from Eq. (7). The solutions 
approach x oo and y I provided that p does not 
pass the singularity at p = po before reaching p ^ pc- 

When pq > pc the system reaches p = Pc with a finite 
time Nc satisfying ^(7 + l)Nc = {r-l + l/xi)''+^ -{r- 
2^7+i_ rjj^j, Hubble parameter vanishes at this point, 
since = Kp/3x from Eq. (6). From Eqs. (9) and (10) 
the differential equations for y and H arc given by 

dy ^ B ^ dH ^ B{l-2y) 

di {r-y)^ ' dt 2(r-?/)^ ' ^ ' 

When ?/ = 1 one has dy/di = and dH/di < 0. Then 
the Hubble parameter becomes negative, which is ac- 
companied by the decrease of y. From Eq. (17) we find 
dH/dt > for y < 1/2, during which H increases. In the 
type I case y and H asymptotically approach zero with 
time-dependence y oc and H cx —t~^. The type II 
case is different because of a time-dependent equation of 
state. In fact when H — and y = we find dy/di = 
and dH/dt > 0. This behavior is clearly seen in Fig. 2. 
Both H and y increase after the system passes the point 
H = and y = 0, which is followed by the maximum 
value of H aX y = 1/2. After that the evolution of the 
universe mimics the previous one, namely H oscillates 
between —1/2 and 1/2 together with the oscillation of 
y between and 1. Hence the universe repeats the cy- 
cle of expansion and contraction without reaching any 
singularities (see Fig. 2). 

When Pq < Pc, independently of p > po or p < po, the 
solutions reach the sudden future singularity at p = po 
in a finite time. 



In standard Einstein gravity the type II singularity 
appears when the pressure density p diverges as p ap- 
proaches some constant value po. For example, this is 
realized when p is given by [5] 



B 



P 



(po - p)'^ 



(14) 



where B, po and 7 are positive constants. This singular- 
ity appears at a finite time as p approaches pp. 

Let us consider the cosmological dynamics in the pres- 
ence of the loop correction. The equation of state is now 
dependent on p, i.e., w = —1 — B/ p{pQ — p)'' . Substitut- 
ing this expression for Eq. (8) by using the relation (7), 
we get 



dx 
dN 



B 



(r - 1 + l/a;)T 



(15) 



C. Type III singularity 

The type III singularity appears for the model 

p^-p-Cp^, /?>1, (18) 

where C is a positive constant. Integrating Eq. (4) for 
the equation of state (18), we find that the scale factor 
is given by 



a = ao exp 



3C(l-/3) 



(19) 



where oq is a constant. In Einstein gravity one has p 
00 and \p\ — !■ 00 in a finite time, but a is finite when 
(3 > 1. Hence this is different from the Big Rip singularity 
at which scale factor diverges. 
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FIG. 2: Evolution of the Hubble parameter and the variable 
y = p/ pc for the model (14) with B = 1, r = po/pc = 2 
and 7 = 2. We choose initial conditions tji — 0.01 and Hi = 
[j/i(l — yi)Y^^- The Hubble parameter H oscillates between 
— 1/2 and 1/2 without reaching the singularity at p — po. 



In LQC Eq. (8) gives 
da; 



1 - 



1 



(20) 



where C = 3Cp^ ^. This is integrated as 

/J-l / \ 



1 



1 - l/x 



1 - l/Xi 



C{f3-1)N. (21) 



Then we get x — *■ oo, y ^ 1 and H ^ a,s N Nc, 
where is given hy C{(3-l)Nc = {l/{l-l/x^)f-^ -1. 
The differential equations (9) and (10) are 



(22) 



Hence one has dy/di = and dH/di < for y = 1 and 
H = 0, which is followed by the decrease of y and H 
(< 0). The evolution of the system is similar to what 
we discussed in the type I case. After the Hubble rate 
reaches a minimum at y = 1/2, H and y asymptoti- 
cally approach H ^ y ~ 0. When y <C 1, in fact, we 
have d^jz/dP « {3C^/2)y^'^ from Eq. (12), which gives 
y oc ^-2/(2/3-1) and H cx -^-1/(2/3-1), Hence the final 
attractor is a contracting universe with p ^ 0, p ^ 
and H —0 as t oo. 



turn cosmology. Non-perturbative quantum effects give 
rise to a correction whose effect depends upon the 
ratio p/pc where pc is of order of Planck density. Typ- 
ically this type of correction is thought to be important 
only in early universe whose energy density is close to 
Pc but it can be also important in future universe if 
(phantom) dark energy is present as observations suggest. 
Note that the modifications we studied are different from 
those given which emerge from the regularization of in- 
verse scale factor operator in LQC and can be important 
below a critical scale factor a* (see for e.g., Ref. [26]). 
These corrections are negligible for a ^ a* and are not 
considered in the present work which deals with late time 
expansion dynamics in LQC. 

There are several types of future singularities which 
appear in standard Einstein gravity. In the type I case 
where p, p and a diverge in a finite time and in the type 
HI case where p and p are infinite but a is finite in a 
finite time, we find that the loop quantum modifications 
generically remove these singularities. The universe tran- 
sits from an expanding branch to a contracting branch 
after the energy density approaches critical value pc. Af- 
ter the Hubble parameter reaches a negative minimum 
when p/pc = 1/2, it increases toward a stable fixed point 
7J = in an infinite time (see Fig. 1). The fate of the 
universe thus dramatically changes on considering loop 
quantum modifications in the standard Fricdmann dy- 



In the type II case where p diverges but p (po) and a are 
finite in a finite time, sudden singularity is not removed 
when Po is smaller than pc- When po > pc, however, 
the Hubble parameter H exhibits an oscillation around 
H ~ (sec Fig. 2). This corresponds to an oscillating 
universe without any singularities. 

We have thus shown that in most cases the future sin- 
gularities are avoided because of the presence of loop 
quantum corrections. Our analysis of the resolution of 
singularities clearly reflects the important role played by 
non-perturbative quantum gravity modifications in order 
to fully understand the dynamics of universe around the 
Planck energy. 
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